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Abstract. We observed the 12 CO(3->2) emission of the emission-line regions Hubble I, Hubble V, Hubble X, Holmberg 18, 
and the stellar emission-line object S28 in NGC6822 with the ESO Atacama Pathfinder Experiment (APEX) 12m telescope as 
part of its science verification. The very low system temperature of 130 - 180 K enabled us to achieve detections in 4 single 
pointings and in a high spatial resolution 70" x 70" map of Hubble V. We compare the spectra with Hi observations, obtained 
with the Australia Telescope Compact Array, of the same regions. In combination with previous multi-line CO observations, we 
perform a preliminary investigation of the physical conditions in Hubble V using a simple LTE model. We estimate the mass 
of the Hubble V region and the H2//co<3^2) conversion factor. Also, we show that Hubble V is located very near the line- width 
versus size relation traced by the Milky Way and LMC molecular clouds. 
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1. Introduction 

We present 12 CO(3— >2) observations of the emission-line re- 
gions Hubble I, Hubble V, Hubbl e X, Holmberg 18, and the 
stellar emission-line object S28 jKillen & Dufourl Il982l) in 
NGC6822 with the ESO Atacama Pathfinder Experiment 
(APEX) 12m telescope 1 . NGC6822 is a Local Group dwarf ir- 
regular galaxy, of type I B(s)m. Its study started with the land- 
mark paper of Hubble (1925), who identified it as a galaxy 
in its own right, external to the Milky Way. The most re- 
cent estimate of its distance, based on ob servations of 116 
Ceph ei d variables, places it at 466 + 20 kp c dPietrzvhski et all 
12004 . IPeimbert. Peimbert & Ruizl d2005l) find a metallicity 
12 + log(0/H) = 8.37 ± 0.09 for Hubble V and 12 + log(0/H) = 
8.19 ± 0.16 for Hubble X, corresponding to roughly half the 
solar metallicity. Star formation proceeded at an almost con- 
stant rate up to the present, except for the central bar region, 
where star-for mation i ncreas ed by a factor of 3 - 4 during the 
last 600 Myr JWvderl l200lh . Its proximity allows us to study 
the different components and phases of its interstellar medium 
on scales of order 10-100 parsec. 
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The detection of the compact molecular clouds associ- 
ated with Hubble V was first reported by IWilsonI {l994). 
Later on, the emission-line regions Hubble I, Hubble V, and 
Holmberg 18, and the stellar emiss ion-line obje ct S28 were ob- 
served in 12 CO(l— >0) emission hv lTsraell (fl997). Moreover, for 
Hubble V, the brightest HII region in NGC6822, detections of 
12 CO(2^1), 1 2 CO(3^>2) 12 CO (4^3), and 13 CO(1^0) have 
been reported Jlsrael et alll2003l I2 CO(1->0) and 12 CO(2->l) 
observations center ed on Hubble X, on the other hand, did not 
yield a detection ( Israe l et all 120031) . These results will serve 
as a comparison for the APEX data presented here. 

In this Letter, we combine the 12 CO(3— »2) line intensi- 
ties measured with APEX with line intensities of other 12 CO 
and 13 CO transitions, taken from the literature, to constrain 
the physical conditions of the molecular interstellar medium 
of NGC6822 using simple LTE models. We also investigate 
the spatial distribution of the 12 CO(3— >2) emission and how it 
correlates with previous high resolution HI observations. 

2. Observations and data reduction 

We observed the 12 CO(3— >2) line towards the star forming 
regions Hubble I, V, X f and Holmberg 18, an d the stellar 
emission-fine object S28 fallen & Dufoutll 19821) in NGC6822 
with the ESO Atacama Pathfinder Experiment (APEX) 12m 
telescope on the nights of 17, 18, 20, 21, 25 and 26 August 
2005 as part of its science verification. We used a Fast Fourier 
Transform Spectrometer (FFTS) backend built by MPIfR 
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Fig. 1. CO(3— >2) spectra of the star-forming regions Hubblel, HubbleX, and Holmberg 18 and the stellar emissio n-line object S28 
in NG C6822, overplotted with the best fitting Gaussian. The black curve indicates the Hi emission as derived by d e Blok & W alter 
(2006). Evidently, the molecular gas associated with Hubble I and Holmberg 18 has the same velocity as the neutral gas. The 
velocity of the molecular gas associated with the stellar object S28 differs from that of the Hi. Hubble X was not detected. 



jKlein et allEo 06). The frequency was centered on CO(3— >2) 
(345.79 GHz) and corrected for the systemic velocity of 
NGC6822. Th e beamsize of the tele scope at this frequency is 
18.2" FWHM jGuesten et alll2006h . We used a 1 GHz band- 
width divided over 8192 channels, corresponding to a veloc- 
ity resolution of 0.1 km s _1 per channel. All targets were ob- 
served with a single pointing, except for HubbleV, which was 
spatially mapped on a rectangular Nyquist 1' x 1' grid. On- 
source integration times were 31.4 minutes (Hubble I), 9. 1 min 
(single pointing of HubbleV), 1.0-10.4 min per pointing of the 
map of Hubble V, or 110.1 min in total, 3.1 min (HubbleX), 
21.9 min (Holmberg 18), and 22.3 min (S28). The rms of the 
pointing model was about 3". Pointing was regularly checked 
and updated on the nearby pointing source W-Aql. We used 
position switching. The reference position was set to +480" 
in R. A. relative to the center position of each source. Dual side 
band (DSB) system temperatures were 128 K (Hubblel), 152 K 
(single pointing of HubbleV), 135-182 K (map of HubbleV), 
133 K (HubbleX), 15 6 K (Holmberg 18), and 176 K (S28) 
(Risache r et all 120061) . The precipitable water vapour was 0.5- 
1.0 mm, corresponding to a T225 of 0.27-0.29, during the ob- 
servations. Calibration errors amount to ~ 15%. In order to en- 
hance the signal-to-noise ratio without compromising the spec- 
tral resolution, the spectra were rebinned to a 0.8 km s res- 
olution and they have been rectified over the velocity range 
[-100, 0]kms''. 

The data reduction was performed with the standard data 
analysis program GILDAS of the 30m IRAM radio telescope. 
Observations related to the same pointing were first added to- 
gether. Afterwards a polynomial of 4th order was fitted to an 



emission line free region of the spectral baseline and subtracted 
off this baseline. Some spectra still showed a residual double 
sinusoidal variation. One of these ripples arose due to a vibra- 
tion of the gore-tex membrane that covers the entrance window 
to the Cassegrain cabin (L. Nyman, priv. comm.) and the other 
due to the vibration of the cold head of the closed-cycle cooling 
machine, which affected the L O coupling of the rece iver, caus- 
ing the receiver gain to vary ( Risache r et all 1200 6). We fitted 
a combination of two sine functions to those spectral baselines 
that contained this sinusoidal variation, omitting the spectral re- 
gion around the CO(3— >2) emission line, and subtracted this off 
the spectrum. The antenna temperatures, T* A , were converted to 
to main-beam brightness temperatures (T m b = T^/r] m b), using 
the main-beam efficiency T] m b = 0.7. 



3. Discussion 

We fitted Gaussians to the detected emission lines in order to 
estimate the peak intensity, T m \, (K), and the integrated inten- 
sity, Iqo = fTmbiv) dv (Kkm s~'), of the 12 CO(3— >2) emission 
lines of the observed star-forming regions (see Table We 
used the best fitting Gaussian and the lcr noise on the spec- 
trum, estimated from the spectral region between -100 and 
km/s and excluding the emission line, to generate 1000 new 
noisy spectra. These were analysed the same way as the orig- 
inal spectrum, allowing us to estimate the lcr errors on these 
quantities. For the non-detected star-forming region HubbleX, 
we give a 3<x upper limit over a velocity width of 6 km/s. 
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Table 1. CO(3— >2) properties of the targeted regions in NGC6822: the peak main-beam temperature, T m b, the velocity of the line 
with respect to the Local Standard of Rest, the line FWHM, and the integrated line intensity, Iqo . 



name 


RA 


DEC 


T mb (K) 


LSR velocity (km/s) 


FWHM (km/s) 


Ico (Kkms- 1 ) 


Hubblel 


19 44 31.64 


-14 42 01.2 


0.07 ± 0.02 


-66.3 ± 0.4 


6.1 ± 1.1 


0.49 ±0.11 


HubbleV 


19 44 52.80 


-1443 11.0 


0.61 ±0.02 


-41.3 ±0.1 


6.0 ±0.2 


3.89 ±0.14 


HubbleV 


19 44 52.80 


-1443 11.0 


deconvolved 


-41.3 ±0.4 


6.0 ±0.2 


6.65 ± 0.59 


HubbleX 


19 45 05.20 


-14 43 13.0 








< 0.30 


Holmberg 18 


19 44 48.93 


-14 52 38.0 


0.11 + 0.02 


-43.6 ± 0.2 


2.5 ±0.4 


0.29 ± 0.07 


KD82.S28 


19 44 57.79 


-14 47 51.5 


0.12 ±0.02 


-57.5 ±0.3 


5.2 ±0.7 


0.67 ±0.11 
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Fig. 2. 12 CO(3 -> 2) map of HubbleV. Each panel shows the 
brightness temperature T m \, as a function of velocity with re- 
spect to the Local Standard of Rest (LSR). Some spectra still 
show some residual variations even after subtracting off a dou- 
ble sinusoidal baseline. Nearest neighbor panels overlap by 8", 
i.e. about half of the APEX beam width at this frequency, so 
that next nearest neighbor panels are roughly independent. The 
CO source associated with HubbleV is clearly resolved in this 
map and is extended towards the north-east. 



3.1. Physical conditions in HubbleV 

We made a preliminary assessment of the physical conditions 
in the CO cloud associated with HubbleV assuming local ther- 
modynamical equilibrium (LTE). In that case, there is one ex- 
citation temperature, T ex responsible for populating the energy 
levels of the 12 CO and 13 CO isotopomers. This need of course 
not be the case in reality, with the higher-/ lines not being ther- 
malized due to their larger Einstein A coefficients. In the fol- 
lowing, we will assume the isotopic ratio X - 12 CO/ 1 3 CO = 60 
( Lang er & Penziasl Il993t ISavage. Apponi. Ziurv in 
which case the optical depths of I2 CO, denoted by Tu, and 
13 CO, denoted by T13, obey the relation Tyi = Xt\^. The cal- 
culated line intensities are coupled to the observed quantities 
by the beam filling factor fc>, which we assume to be the same 
for the 12 CO and 13 CO emission. Furthermore, we will assume 



Table 2. Comparison of the LTE model with the measured in- 
tensities of the 12 CO and 13 CO transitions. All intensities have 
been corrected for beam dilution using Q S0U rce = 209 arcsec 2 . 



transition 


measurement (K) 


model (K) 


12 CO(1^0) 


1.72 ± 


0.11 


1.73 


12 CO(2->l) 


1.69 ± 


0.08 


1.66 


12 CO(3->2) 


1.56 ± 


0.05 


1.57 


12 CO(4->3) 


1.50 ± 


0.09 


1.48 


I3 CO(1^0) 


0.11 ± 


0.07 


0.11 



that both the 12 CO and the 13 CO emission arises in the same 
region so that Q SOU ice, the solid angle spanned on the sky by the 
CO source, is the same for both isotopomers. This is to keep 
this preliminary modeling as simple as possible since there is 
no physical reason why fi, and f2 source should be the same for 
all transitions. We can use this source solid angle to correct the 
observed emission line brightness temperatures for beam dilu- 
tion using the relation T' mb = r m b(O source + O k J /Q source , with 
£2beam the beam solid angle. The main-beam brightness temper- 
ature of an observed transition can be written as 



1 



1 



■,hvt/kT a 



1 



e hVilkT mh _ I 



(1) 



with v, the frequency of the transition, t, its optical depth, and 
Tcmb = 2.725 K the background radiation temperature. 

Using a non-linear minimisation routine, we simultane- 
ously fitted a Gaussian model for the spatial distribution of the 
CO emission of HubbleV, convolved with the APEX beam, 
which constrains the source solid angle Q 80uroe and the pa- 
rameters r ex , ft,, and Tyi to the I2 CO(3— >2) map, presented 
in fig. |2 and to the 12 CO(3— >2) brightness temperature mea- 
suredby us, and the I2 CO(1-*0), 12 CO(2^1), 12 C O(4^>3), and 
13 CO( l— >0) brightness temperatures presented in llsrael et alJ 
(2003). We then used the best fitting values for Q SOU i-ce, T ex , fi» 
and tu to generate 10000 mock data sets with added Gaussian 
noise on each of the observed quantities, using the measured 
1 <x uncertainties on the measured quantities as estimates for the 
dispersions of each of the noise distributions. These mock data 
sets were analysed the same way as the original set, allowing 
us to estimate the 1 cr errors on the derived quantities. This way, 
we find that the parameter values O s ouree = 209 ± 50 arcsec 2 , 
T ex = 49 + 27 K, Tn = 3.7 ± 2.3, and f b = 0.04 + 0.03 pro- 
vide the best fit to the whole data-set. There is a large degree 
of degeneracy between the parameters of this model, e.g. be- 
tween T ex and fb. This is reflected by the very large errorbars 
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Fig. 3. Spatial distribution of the CO(3— >2) emission of 
Hubble V at 18" resolution (contours) plotted over an Ha im- 
age (greyscale). The square indicates the region mapped by our 
APEX observations. The contour values are 1 (which corre- 
sponds to 3<t), 1.5, 2, 2.5, 3.0, and 3.5 Kkm s . The CO emis- 
sion is not centered on the Ha emission. The posi t ions o f the 
CO clouds MCI and MC2, identified bv Iwilsonl dl994h . are 
indicated in the figure. 



on these quantities. Still, the minimisation routine converges to 
the same solution independent of the starting point of the min- 
imisation which proves that the minimum of the x 2 is well de- 
fined. Moreover, this temperature estimate agrees with the dust 
temperature derived from the ratio of the 60yU and 100/i IRAS 
flux densities, / v (60) = 7.89 Jy and / v (100) = 11.81 Jy, of 
HubbleV, 7d us t ~ 40 K. This estimate was derived assuming a 
single temperature component and a A emissivity law. Given 
the apparently rather high temperature of this CO emission 
cloud, observations of higher-/ transitions, e.g. with FLASH 
or the future SIS heterodyne receivers, are required for a more 
precise assessment of its physical properties using more sophis- 
ticated LVG models, taking into account non-LTE effects. Also, 
some of the publishe d high-./ transition temperatures, such as 
12 CO(4->3) value of llsrael etal] d2003h . may be affected by 
the small area that was mapped. If some of the emission was 
missed, this may lead to an underestimation of the brightness 
temperature of these transitions. 

Using the 12 CO(3->2) FWHM line-width, AV, in km s 1 , 
and the radius of the emission region, R, in parsec, we can also 



estimate the virial mass of Hubble V as M v 



in solar masses (MacLaren, Richardson, Wolfendale, 
For AV = 6.0 km s -1 , R = 
6.1" = 13.9 pc, we find M vil - * 9.5 x 



190R(AV)\ 
1988). 



source /2 . 3 5 5 — 

1O 4 M . Using 

the relation M rimt = 1.27f v (100 )£> 2 (exp(144g/r) - 1)M 
( Bos elli. Leaueux. GavazziL 20021) for the dust mass, with D 
the distance in Mpc, we find Md us t ~ 130 M . Using the 
metallicity-depend ent gas versus du st mass relation (eqn. (6) 
in Boselli, Leaueu x". Gavazzil ( 120021) '). this yields M gas w 9 x 
1O 4 M , in good agreement with the virial mass. This is much 
less than the estimate of the total gas mass M gas = 10 + 5 x 



10 5 M of llsrael et al.l d2003l) . This is most likely because we 
are measuring quantities (radius, velocity dispersion) that per- 
tain only to the CO emission region : any mass distribution out- 
side t his region, which was taken into account bv llsrael et alJ 
il2003t) . would not have a very large effect on the virial mass 
estimate. Hubble V is also very near the line-width versus size 
relat ion traced by the Milky Wa y and LMC molecular clouds 
jHeikkilaL fl99l iMalonevt 1 1 990l) . We estimate the 12 CO(3^2) 
to H2 conversion at 5.4-5.8 x 10 20 crrT 2 (Kkms -1 ) -1 , depend- 
ing on whether the deconvolved or the observed brightness is 
used. 



3.2. Comparison with previous observations 

In Figures [H and 13 w e plot the (unsealed) HI spectra derived 
bv Ide Blo k & Walter (2006), of the same regions as our ob- 
servations on top of the 12 CO(3— >2) spectra. One can see that 
apart from S28, both emission lines are coincidental although 
the HI emission is systematically broader than the 12 CO(3— >2) 
emission. In Figure [3] we plot the spatial distribution of the 
12 CO(3— >2) emission of HubbleV at 18" resolution . We find 
a diffe rent morphology than the one derived by llsrael et all 
(2003), however it is in accordance with their [CII] emission, 
found in the same paper. A possible cause might be the higher 
system temperature of 2460K during their observation, produc- 
ing more noise which might cause a shift in the spatial dis- 
tribution. The main emission peak in our map corr esponds to 
the molecular cloud MC2, first detected bv lWilsonl Jl994l) : the 
eastward extension corresponds only very roughly to MCI. 

With our observations we prove that APEX is very suited 
for deriving spatially extended, high signal-to-noise maps of 
emission-line regions in Local Group dwarf galaxies, where 
one can achieve a spatial resolution of a few tens of parsecs. 
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